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Abstract. For a maximally entangled eigenstate of a system of two non-interacting identical one 
dimensional harmonic oscilators, at the semiclassical level, it is not obviously true that a nonlinear 
interaction with one of the subsystems leaves the reduced semiclassical Wigner function of the other 

CsJ ■ subsystem unaffected. Once stated, we advance some wild speculations regarding this seeming puzzle. 
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! When Einstein, Podolsky and Rosen explored the "strange" behaviour of entangled states, their 

main goal was to set forth an argument concerning the ontology of quantum mechanics. Only quite 
I> ' later, did Bell show that some of the discussion had a definite empirical bearing, namely, on 

■ matters related to the nonlocal character of quantum mechanics. After the empirical confirmation of 



O ' this latter [|1| , emphasis has been placed on the "peaceful coexistence" |]T5| of this kind of nonlocality 



^ I and the partial order of spacetime events derived from the principle of relativity 
^ ■ In terms of the Wigner |]T^ function W, the cartesian Weyl representation of the density 

O ■ operator p as a real function on euclidean phase space, this question is addressed in the following 

^ ! way: let Xi and X2 denote the phase space points of two subsystems whose total system is in an 

^1 entangled eigenstate of the hamiltonian Hi + H2 , with corresponding Wigner function Wq{xi,X2)- 

[ Also, let H2 be a new hamiltonian, whose Weyl symbol H2 generates a flow X2(0) — ^ X2{t), and let Wt 

■ denote the evolution of Wq effected by this new hamiltonian. Further: let Wq{xi) = J Wo{xi, X2)dx2 

■ and W^{xi) = J Wt{xi, X2)dx2 be the corresponding reduced Wigner functions, which encode all the 
>- ! empirical information of subsystem 1 alone, and denote their difference by 5W^{xi,t). 

\ In this context, conservation of the relativistic partial order of events is implied by 

. 9^. ; (1) = , 

which is always true in full quantum mechanics {%= 1), due to its linear structure [0]. In a classical 
setting {h= 0), on the other hand, a sufficient condition for the validity of ([l|) is given by 

(2) 1^0(^1, 2:2(0)) — ^ Wt{Xi,X2{t)) = W^{xi,X2m . 

Years ago. Heller noticed that the generic expression corresponding to @j is not always valid at 
the semiclassical level [h —>■ 0~^), a fact which has recently been reinterpreted in a more geometrical 
way [Q. Accordingly, for a generic semiclassical Wigner function W, (H) is generally true only if if 2 
is quadractic. For a nonlinear interaction, is valid only for points X2 close to the Lagrangian leaf 
corresponding to the classical limit of W, whenever its semiclassical expression W is an oscillatory 
function of X2, as is usually the case for pure states • 



The purpose of this letter is to state the following puzzle: at the semiclassical level, it is not obvious 
whether (P is valid in general. To see this puzzle, we focus on a particularly simple semiclassically 
entangled state, namely the (anti) symmetric state of two identical 1-dimensional harmonic oscilators, 
under a small nonlinear interaction H2. We apply the results of [H] to ask whether, in a suitable limit, 
(0) might not hold. Finally, we emphatically stress the highly speculative nature of this letter. 
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2. THE SETUP 

We study the semiclassical dynamics of an eigenstate oi H = Hi + H2 , where Hi = {pf + gf)/2, 
and of -P12 : 1 ^ 2. This state is: = \em, 1) ® |e„, 2) ± |e„, 1) le^, 2) , where the ±1 refers 

to the eigenvalue of -P12 and where Cm+n = + e„ is the eigenvalue of H , with = (m + l/2)/i , 

= ('^ + l/2)h and where the energy difference is A = — e„ = Zfi, with m — n = I > 0. 

The corresponding Wigner function factors as 2Wm+n{xi,X2) = Wm{xi)Wn{x2) + Wn{xi)Wm{x2) i: 
M^{xi)M^{x2) ± M^{xi)M^{x2) , where is the Moyal fll]] function, or cartesian Weyl repre- 



sentation of the transition operator |em)(e„|, satisfying = M^. We write this function in polar 
{x = {r,6)} coordinares but emphasize that we are not considering the action-angle Weyl represen- 
tation, which is not equivalent 0. In taking the semiclassical limit 0, we formally let ^ ^ 0"*" . This 
is equivalent to letting m, n — > 00 , while keeping Cm, e„ < cxd , or to varying an "effective" h. 

Our analysis is based on the semiclassical expression for the Moyal function which is not corrected 
by uniform approximations on the caustics. Ours is a particularly simple case of the general expression 
carefully treated in [|T^] and is given by (we use curly letters to indicate semiclassical expressions): 

where (j)m,n{f) is half of the symplectic area between the circle with radius a/ (2m + 1)^ centered at 
the origin and the circle with radius {2n + 1)H centered at {2r,6) , and Dm,n{f) = l^m ^ ^n\^^'^i 
where x^ , x~ are the phase space velocity vectors of the intersecting circles. These intersections are 
real as long as r_ < r < r+, where r± = (^ (2m + l)h± ■>/ (2n + l)^)/2 , and the denominator goes 
to zero when the circles become tangent, i.e. r = r± , which are the caustic lines. Outside this ring, 
A^J^ = , in this crude expression. Naturally, the corresponding formula for Wm is obtained from 
(0) with m = n. By performing uniform approximations as the intersection points coalesce, can 
be expressed in terms of Airy functions which do not blow up and decay exponentially outside the 
ring. Still using (|^), the sum of the cross terms in Wm+n is expressed as 

(A\ nlU^f U = - ^2))cOg(0m,n(n)/^ - 7r/4)cOg(0^,„(r2)//l - TV / 4) 

KeiJvi^^[Xi)jvi^,[X2)i {n^h/2)D^,^{n)D^,^ir2) 
Of course, / Re{M"^{xi)M';^{x2)}dx2 = and so = (Wm + W„)/2 , as it should be. 

3. THE PUZZLE 

The seeming puzzle to be laid forth below follows from the fact that i^ oscillatory function, 
that is, (0) has an essential singularity at ^ = 0. This is the reason why (|^) fails to hold when ^ — 0"*" 
as can be seen by directly plugging (^) into the Moyal bracket. Its successive derivatives then bring 
negative powers of h which alter the usual ^-expansion obtained for smooth functions. Recently 
T^ , a simple geometrical prescription for a more correct semiclassical propagation of such functions 



was obtained, in agreement with another old result derived by Berry and Balazs according to 
which one should propagate the whole family of curves, or lagrangian submanifolds, associated with 
a quantum state and then re-evaluate the semiclassical expressions accordingly. 



The new result obtained in [14| somehow simplifies this procedure by showing that, equivalently, 
one can obtain a more correct semiclassical propagation of an oscillatory function like A^^(x) by 
classically propagating the tips of the chord centered on x instead of classically propagating the 
argument x itself, as in the Liouville case (|]). Of course, Lioville propagation is still semiclassically 

^There are different phenemenological ways of approaching the classical limit. The one considered here, large 
quantum numbers or small effective Planck's constant, should not be identified with the limit of many degrees of 
freedom, or large complexity, for which most present theories of decoherence apply. 
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correct for smooth functions, when all chords are null, or if the flow is linear, as discussed in more 



detail in |jT^. Since A45^(x) is not smooth, we might expect surprises when under nonlinear interac- 
tions. And these might be related, in some way or another, to the old strange fact that for generic 
pure-state Wigner functions under nonquadractic hamiltonians, if the interaction is first performed 
in the full quantum regime and then the classical limit is taken, the final result generally appears to 
be rather different from first taking the classical limit and then interacting. 

Therefore, let us now suppose that subsystem 2 interacts nonlinearly for an interval of time t. For 
simplicity, we state the puzzle when the interaction is either cubic H2 = or quartic H2 = ttg|/4 

and denote by e = at the strength of this interaction. After rewriting (^) as a sum of oscillatory 
functions of 4>mni^2,02) = 4'm,nij'2) ± , it is uot too difficult to see that, within the stationary 
phase regime {h —>■ 0"*"), the effect of a small interaction is, in a first approximation, the equivalent 
to having only phase shifts 0^,„(r2, 62) <Pti,ni.'r2, O2) + 5(/>^,„(r2, 6*2). 

More precisely, to get these phase shifts we first apply equation (16) in [jl4| to the chord ^^ .^ = 
— J(90^„(x2)/5x2), centered on X2, and get the "midpoint phase difference" 0m,n(^2) ~ </'m,n(^2), 
where x\ midpoint of the hamiltonian flow of the tips of the chord ^. Then we note that, for 

small e, the total phase difference acquired for the hamiltonian flow of X2 is approximately half of 
'Pm,ni^2) ~ 4>m,ni^2)- The amplitude being approximately covariant under the hamiltonian flow, for 
small e, this determines the looked-for phase shifts as 

50i„(r2,02) ^ (2e/3){szn(^2)[e^+n - (A/2r2)2 - r^jV' ± 00^(^2) [A/2r2]}3 , 
for the cubic case, while for the quartic case we get 

50±,Jr2, 62) ^ {2e)r2Cos{e2){stni92)[e^+n - (A/2r2)2 - 4]'/^ ± 003(^2) [A/2r2]}3 . 

The next step is to compute / Re{J^^{xi)A4''^ {x2)}dx2 and check if it vanishes, as before. But 
note that the phase shifts have broken the symmetry of Re{Ai^{xi)Ai'!^ {X2)} with respect to 62 
and it is not obvious now whether the integral vanishes or not. In fact, we must now perform the 
new integration f Re{A4^{xi)J^^ {x2)}dx2 by stationary phase and it is not difficult to see that 
the critical points of 0^^(r2,6'2) + 50^ „(r2,6'2) are generally complex. This means that we must 
generally perform a 2-dimensional steepest descent evaluation of the integral, carefully deforming the 
whole ring into the complex 2-plane, counting all the significant contributions, adding them all up... 
This is a rather delicate computation, but it is not unreasonable to suppose that, instead of zero, we 
should at least get an exponentially small result Q], as in a "tunneling" effect 0. 

This being an asymptotic analysis, however, we still need to make sure whether the result of this 
stationary phase computation implies that SW^ is really different from zero, to lowest order of h 
(remember that (^ is only the first term of an ^-expansion), in some region of the positive (A, e, h)- 
space close to the h = plane. Also in this region, we must analyse the behaviour of / yVm{x2)dx2 and 
J Wn{x2)dx2 under similar phase shifts to find out how unitary such a propagation of actually is 
(note that J J Re{J^'^{xi)J^^{x2)}dxidx2 vanishes). Then we must compare all the above results 
and further analyse and interpret them (for instance, in possible relation to Stokes phenomenon). 

Finally, to properly account for stronger interactions (e not so small), closer resonance (A small), 
or to actually get more precise estimates for (5W^ (any e or A), we must leave the crude expression 
(0) and proceed much more carefully with the uniform approximations (this is the really difficult 
part), since non-Liouville propagation of the amplitude is so much more relevant near caustics. 

It is very important to emphasize that any possible nontrivial answer for 5VV^ would only be an 
asymptotic one. Such a possibility seems due to the 6'2-symmetry breaking of the whole semiclassical 

^ For a much simpler analogy, consider the integral X = J^^ cos{l0)d9 . Writing 16 = (l>{6)/h , where (/){6) ~ X6, 
then under the phase shift (f> ^ (j) = cf) + 6(f), where Scf) = e 9^/3 , we have that I ^ X is asymptotically (ft — > 0+) 
mapped into the Airy integral, that is, X « T:{h/ey/'^Ai{y) , where y = Ae^^/'^ft^^/"^. 
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Wigner function yVm+n{xi,X2), effected by the nonlinear interaction with H2 , and the consequences 
of this symmetry breaking should become more explicit as the amplitude corrections are included. 
But, actually, this is precisely the very hard crucial point that needs to be fully investigated. In 
other words, it is quite likely that only a careful (and very difficult) analysis based on more accurate 
semiclassical approximations can produce clear answers, even because the accuracy of the chord 
expression (^, from which we started, is in fact worse than exponential. 

Then, z/such a possibility is indeed correct, we could conjecture there would be an optimal way of 
approaching the classical limit, in the positive (A, e, ^)-space, leading to a maximal 6W^. Actually, 
combined with the singular nature of the semiclassical limit, the presence of so many independent 
parameters is what seems to open this problem in such a manner. And for more degrees of freedom, 
issues from the classical (non)integr ability of a subsystem, with respect to the nonlinear interaction, 
might well play a further significant role. However, it is clear that nothing in the previous discussion 
proves or even suggests that a nontrivial (5W^ can probably happen. In this letter, we have merely 
pointed out to an apparent possibility that this somehow might come to pass. 

4. SPECULATIONS 

Final words on this puzzle from the phenomenological, empirical and philosophical points of view. 

Phenomenologically, how could it ever be possible for (|l|) to fail at all, since it is valid in both the 
classical and the full quantum regimes ? Well, we could argue that (|I]) is valid in both regimes 
for different reasons. In the classical case, translates the fact that no quantum correlations have 
"survived". But this can be seen as an effect of the correlations becoming "random", which is best 
expressed by the oscillations in the Wigner function becoming so high that they average to zero 
everywhere but on the classical leaf. In the full quantum case, by contrast, the correlations are fully 
structured, in other words, they have a definite "form", which is usually expressed by the concept 
of a definite normalized vector in the Hilbert space. Such a vector can only be preserved by linear 
transformations, from which follows. Therefore, the puzzle ammounts to a question of whether, 
at the semiclassical level, the linear structure of quantum mechanics could feel some disturbing effects 
of a nonlinear classical dynamics and, if so, how they could show up 0. 

Empirically, if a suitable limit where 7^ could ever be found, would this mean the possibil- 
ity of faster-than-light-communication signals? Here, we should stress that these sub-hypothetical 
"/a/zco" signals would be fundamentally threshold phenomena and so a deeper understanding of this 
threshold between the quantum and the classical regimes could be required. In other words, the very 
existence of stable semiclassically entangled states would be at question. On the other hand, even if 
granted their existence for some special systems, in a suitable limit of nontrivial 5W^, it could still 
be a formidable, pehaps nearly impossible task to actually measure such signals, even though some 
direct measurements of Wigner functions (if ever needed) have already been claimed |]10| . 



Philosophically, what would be at stake here is, among other things, the principle of relativity. 
But in this respect, it seems interesting to point out the usually overlooked fact that, as originally 
formulated ||^, the relativistic partial order presupposes events which are correlated by an exchange 
of energy. In other words, which are correlated by "material" causes. Such is not the case for the 
quantum correlations treated here, which could pehaps be described as "formal" correlations. Maybe, 
these two kinds of correlations being essentially distinct, their principles need not be the same. 



■^Examples of phenomena that appear in one or the other but not in both regimes are easier to come by. On 
the other hand, the non-geometrical phenomenon of the glory fl^ can be seen as a phenomenon that vanishes both 
"classically" (ray regime) and "fully quantum" (long wave regime) but not "semiclassically" (short wave asymptotics) . 

^In a certain way, the two important factors for a Bell's telephone, namely the distinction between pure and mixed 
states and between linear and nonlinear evolution, show up in the semiclassical propagation of Wigner functions. 
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